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ABSTRACT 

We present a new upper limit to the 21 cm power spectrum during the Epoch of 
Reionization (EoR) which constrains reionization models with an unheated IGM. The 
GMRT-EoR experiment is an ongoing effort to make a statistical detection of the power 
spectrum of 21 cm neutral hydrogen emission at redshift z ^ 9. Data from this redshift 
constrain models of the EoR, the end of the Dark Ages arising from the formation of 
the first bright UV sources, probably stars or mini-quasars. We present results from 
approximately 50 hours of observations at the Giant Metrewave Radio Telescope in 
India from December 2007. We describe radio frequency interference (RFI) localisation 
schemes which allow bright sources on the ground to be identified and physically 
removed in addition to automated flagging. Singular-value decomposition is used to 
remove remaining broadband RFI by identifying ground sources with large eigenvalues. 
Foregrounds are modelled using a piecewise linear filter and the power spectrum is 
measured using cross-correlations of foreground subtracted images. 
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1 INTRODUCTION 

Between the recombination epoch at 2: ~ 1100 and the 
first round of star formation at z ~ 10 the Universe was 
filled with neutral hydrogen. This neutral gas is thought to 
have produced 21 cm hyperfine emission with an effective 
continuum br ightness temperature between -500 and 30 mK 
iFurlanetto et al. 2006) using z = 8.6 a.nd W MAP 7-year 
cosmological parameters (|Komatsu et al.|[2010l ). which have 
formal error bars of roughly 10 per cent. As the first stars 
formed, beginning at local peaks in the matter density, the 
hydrogen gas was locally ionized. This is the start of a period 
of cosmic evolution called the Epoch (or Era) of Reionization 
(EoR). Over time the ionized cells grew and overlapped, cre- 
ating a patchwork of ionized and neutral cells. This general 
patchy topology is well motivated by theory and simulations 
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though the exact properties are poorly constrained. Even- 
tually, by z ~ 6, the ionization of the broadly distributed 
material was complete, leaving only rare pockets of neutral 
gas. 

In addition to directly constraining the redshift of 
last scattering zis ~ 1100, Cosmic Microwave Background 
(CMB) polarization data have been used to constrain 
the redsh ift to the range 8.0 < z < 12.8 (WMAP7; 
[Komatsu et al.ll2010l l under the assumption it was instan- 
taneous, although evidence from the CMB indicates tha t 
reionization was an extended process (|Dunklev et al.ll2009l '). 
There is also substantial information on the neutral content 
of the intergalactic medium (IGM) since z ~ 6 via observa - 
tions of the Gunn -Peterson trough |Gunn fc PetersonI 19651) 
in quasar spectr a ( Diorgovski et al.ll200ll : lBecker et al I I2OOII : 
iFan et all I2OO6I : IWillott et all |2007| V Lyman-alpha (Ly-a) 
absorption provides a sensitive probe of neutral gas density 
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and can be used to constrain reionization (e.g., iFan et al] 
|2002| ). With these two techniques we can observe the start 
and the end of the neutral era. However, neither CMB nor 
quasar observations allow detailed examination of the reion- 
ization era itself. For this it has long been proposed to use 
21 cm fluctuati ons in the brightness temp e rature of neu- 
tral hydrogen jSunvaev fc ZeldovichI Il975l : iHogan fc ReesI 
1 19791 ). The 21 cm power spectrum, resulting from a com- 
bination of the patchy ionized and neutral medium and 
the underlying mass power spectrum, is g enerally con- 
sidered one of the mo st promising signals (|Scott fc Reed 
Il990l : IZaldarriaga et al. 2004) and much attention has been 
paid in the literatu re toward designing suitable experiments 



e toward designing suit a Die experiments 
Morales fc Hewitt) |2004| : iMoraled [ioosi : 



to detect it (e-r^ . .- - - 

iBowman et~al]|2006l : iHarker et al.ll2010l V 

Here the observing frequency is 1420 MHz (1 -I- z)~^ , 
in the very high frequency (VHP) range of the radio 
spectrum near 150 MHz. Several programs are underway 
to study the EoR i n the VHP band inclu ding LOPAbQ 
l|Kassim et all |2004|: iRottgering et all l2006l: for the EoR 
case see e.g.. 'Zaroubi fc Sillj bOOalHarker et all HoiO), 
M WAPI (ILonsdale et al.. ,2009; for the EoR case see, 



[Morales et al.'ll2006l: iBowman et al.ll2006l : iLidz et al.|[200;_. 
PAPEBl^l (IParsons et al.l I2OI0I ). 21CMAFI fa.k.a. PaST: 
iPeterson et al.ll2004l ), in addition to the current GMRTrl 
program, and the first goal of such efforts is to determine 
the redshift at which roughly half the volume of the uni- 
verse is ionized. Because the 21 cm emission is an isolated 
single line, these observations can provi de three dimensional 
information over a large reds hift range (|Loeb fc Zaldarriagal 
I2OO4I : iFurlanetto et al.l I2OO6I I. allowing a broad search for 
the half-ionization redshift. In the future 21 cm tomography 
has the potential to add strong constra i nts on cosmological 



parameters (e.g . , iMcQuinn et al 
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20061: ICoorav et al.ll2008l : 
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The brightness temperature of the 21 cm line relative 
to the CMB is determined by the ionization fraction of 
hydrogen and the spin temperature of the neutral popu- 
lation, which is in turn governed by the background radia- 
tion a nd the kinetic temperature of the gas (Purcoll & Field 
Il956l : iFieldl I1959I : IFurlanetto et aLll2006l '). Reionization re- 
quires a minimum expenditure of 13.6 eV of energy per hy- 
drogen atom. In contrast, Ly-a photons, when absorbed by 
a neutral atom, are quickly re-emitted. Ly-a photons are 
said to undergo 'resonant scattering' and each (multiply- 
scattered) photon can affect many atoms before cosm i c red- 
shifting makes them ineffectiv e (IWouthuvsenI Il952l : iFieldl 
ll959l : IChuzhov fc Shapir"oll2006h . This means Ly-a pumping 
of the hyperfine transition requires only about 1 per cent 
of the UV flux required for ionization. Assuming a gradual 
increase in UV flux with time, well before flux levels for 
reionization are reached, Ly-a pumping will couple the spin 
temperature to the kinetic temperature of the gas. If there 
was no source of heat at that era other than a weak UV flux, 
the gas kinetic temperature must have been at its adiabatic 
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expansion value 1.7K at z = 8.5, and neutral gas will pro- 
duce a signal due to absorption of CMB photons in excess of 
stimulated emission and ionized str uctures would be seen as 
low-brightness regions on the sky (|Chen fc Miralda-Escudj 
2004). In such a cold-gas model the brightness temperature 
of the neutral gas against the CMB can be as low as -500 mK. 

A small fraction of the mass will have collapsed into 
minihaloes, which h ave a tem perature higher than this adi- 
abatic temperature jShapiro et al. 2006.'). This fraction may 
even be s ubstantially smaller due to non - perturbative veloc - 
ity flows (|Tseliakhovich fc Hiratal I2OI0I : loalal et al.l I2OIOI ') . 
No study was found on the impact of these non-linear effects 
on the Ly-a pumped IGM temperature. At linear order, the 
positive and negative over and under densities cancel, and 
the non-linear collapse fraction is still small, so we expect 
the realistic value to be similar to the adiabatic prediction. 

Alternatively, X-rays from supernovae or QSOs 
might have heated the IGM above the CMB temper- 
ature of pa 30 K before reion ization ijMadau et al.l Il997l : 
IChen fc Miralda-Escudil2004f ). X-ray heating between 30 K 
and 10 000 K will result in a largely neutral but warm IGM, 
which would be seen in emission. In the limit that Ts 2> 30 K, 
the volume emissivity becomes independent of tempera- 
ture. Patchy X-ra y heating can also result in large angu- 
lar scale structure l|Alvarez et al.ll2010l '). The sky brightness 
temperat ure of the neut ral gas has an asymptote at ~ 30 mK 
((Furlanetto et al.lliooel '). 

The cosmic luminosity of X-rays at 2 : ~ 9 is not 
known and difflcult to estimate (see e.g., iDiikstra et al.l 
|2004 ISalvaterra et al. l l2007l : lGuo et"Zl l200gl. If the rate of 
core collapse SNe at high redshift matches that of today, 
the X-ray output at 2 ~ 9 would have been sufficient to 
raise the IGM temperature above the CMB temperature at 
the onset of reionization. However, the mechanism of and 
factors affecting core-collapse are not known and most nu- 
merical models app ear to generically not result in SNe at all 
jMezzacappal [2005[ ) . It is possible that at high z, core col- 
lapse SNe were not as abundant as today, and the IGM was 
still in absorption during the EoR. We therefore consider two 
limits: one where the IGM is still cold and Tb = — 500mK, 
and another where it is heated above CMB and Tb = 30 mK. 
A general p arametrization of these s cenarios was recently 
proposed bv lPritchard fc Loebl (|2010l ). 

The GMRT-EoR p roject uses the Giant Me trewave Ra- 
dio Telescope (GMRT; I AnanthakrishnanI 1 19951 ) near Pune, 
India, to make measurements of the power spectrum of the 
neutral hydrogen signal with the hope of characterizing the 
structure in the range 8.1 < z < 9.2. The FWHM of the 
GMRT primary beam at 150 MHz is 3.3° which provides a 
cylindrical comoving survey volume of (280 /i~^Mpc)^, with 
about equal dimensions in three directions. The primary sen- 
sitivity comes from the compact central core which is con- 
tained within ~ 1 km, or 30 dish diametres, which gives im- 
ages with ~ 30^^ resolution elements. As recorded, the data 
have high spectral resolution along the line of sight, but we 
bin the data for comparable transverse and radial resolution. 
Each data cube has ~ 30^ ~ 30000 resolution elements so 
even though the signal to noise ratio of each element is less 
than one, the spatial power can be measured to 1 per cent. A 
statistical measurement of the power spectrum provides the 
best hope of pinning down the EoR half-ionization redshift. 

In section[2]we describe the data and analysis, including 
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RFI removal and foreground filters used. In section |3] we ex- 
tend the analysis to cross-correlations of multiple nights, and 
a measurement of the power spectrum. We then conclude in 
section |4l All distances are in comoving coordinates. 



2 OBSERVATIONS 

2.1 GMRT and Data Description 

The GMRT is a radio interferometer consisting of 30 an- 
tennas, each with a diametre of 45 m. Fourteen of these are 
arranged in a dense central core within 1 km which allows 
the high brightn ess sensitivity re quired to search for the 
dim EoR signal (|Pen et al.l |2009| ) . The longest separation 
between antennas is about 25 km. For this experiment the 
telescope was operated at 150 MHz with a 16.7 MHz band- 
width. For 21 cm (1421 MHz) emission, this frequency range 
probes a redshift range of 8.1 < 2 < 9.2. The highest angu- 
lar resolution at this observing frequency at GMRT is about 
20 arcseconds. 

Each antenna provides a pair of left and right circularly 
polarized signals, which are passed through a new signal 
processing system, built in part for this project , called the 
GMRT Software Backend rGSB: lRov et al.ll201(]| '). The final 
recorded visibilities have a resolution of 7.8 kHz and 1/4 s in 
frequency and time, respectively. 

The field is centred on the pulsar BQ823+26. This pul- 
sar has a period of 0.53 seconds and an average flux of 
350 mjy at 150 MHz (Hobbs et al. 2004). It is situated in 
a relatively cold part of the sky at galactic latitude ~ 30° 
with few nearby bright sources. The on-pulse flux is about 
6 Jy, brighter than all other sources in the field, making it a 
g ood calibra t or. T he calibration of this data was described 
in lPen et al.l (|2009|) and is summarized below. The positions 
of the brightest sources in the field are shown in Fig. [1] 

The data are folded in time into 16 'gates' such that 
the pulse from the pulsar is contained within one gate. The 
period of the pulsar is much shorter than the time-scale at 
which ionospheric fluctuations dominate, so by comparing 
the 'on' gate with the neighbouring 'off' gates, everything 
in the field that is constant over the period of the pulsar can 
be removed. This includes sky sources and most radio fre- 
quency interference (RFI). By comparing the pulsar signal 
received by each antenna, the relative system gain of each 
antenna can be calibrated. This technique allows calibration 
of both phase and polarization in real time, with 100 per 
cent observing efficiency, since no additional time is needed 
for phase or flux calibration. Errors in clock synchronization 
sometimes cause the pulsar pulse to straddle two gates. In 
these cases it is necessary to include the flanking gates when 
identifying the pulsar signal. This is not ideal since the pul- 
sar signal is diluted over multiple gates, but the correction 
is limited to only the portion of observations that require it. 

Since the pulsar amplitude varies from pulse to pulse, 
the absolute system gain needs to be measured separately. 
This is done using a noise injection system, which the GMRT 
software backend decodes to calibrate the absolute gain of 
the system, while a sky radio source is used to transfer the 
noise source calibration on to the sky. The primary fiux 
calibrator used was 5C 7.245, a radio galaxy at z ~ 1.6 
l|Willott et al.ll200ll ). located within the field of view (see 
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Figure 1. The 12 brightest sources in the field used for these 
observations. Pulsar PSR B0823-I-26 is at the centre of the field. 
The circle denotes the half-power diametre of the main beam. 



Fig. [1} . However, since this radio galaxy is an extended ob- 
ject, it can only be used as a calibrator for short baselines 
where its structure is not resolved. For antennas in longer 
baselines we determine the relative gain calibration using 
the pulsar. 



2.2 Data Analysis 

Data from 2007 December 10, 11, 14, 16, 17, and 18 are 
included in the current analysis. A software pipeline has 
been developed to automate the calibration and interfer- 
ence removal steps described below. The current configura- 
tion takes approximately 11 hours to process one hour of 
data on the CITA Sunnyvale computing cluster, or 20 hours 
on the National Centre for Radio Astronomy (NCRA) HP 
computing cluster in Pune, India. The NCRA facility can 
process four hours in parallel, completing one night of ob- 
servations in slightly less than two days. While the Sunny- 
vale cluster has the capacity to process as much as 12 hours 
of data at once, storage capacity, shared usage, and other 
bottlenecks reduce this significantly. 

One of the limiting factors for measurement of the 
EoR signal is broadband radio frequency interference (RFI), 
which dominates the signal at 150 MHz. In addition to the 
standard procedure of flagging bad data, we also use a sin- 
gular value decomposition (SVD) to remove broadband RFI 
from the data, and also to identify and remove interference 
sources. We believe this approach is unique among other 
RFI mitigation strategies in the radio community. 

Narrow line interference was removed by masking points 
in each frequency bin with an intensity above some thresh- 
old. This was done twice in the data reduction pipeline. In- 
put data are initially flagged with a threshold of 8(7 on a 
Gaussian scale, and then again at 3(t after broadband in- 



4 G. Paciga et al. 











1 





















Figure 2. Representative stages in tlie data reduction pipeline. 
In eacli panel, the horizontal axis is frequency, covering approx- 
imately 1 MHz, and the vertical is time, increasing downward, 
and covering approximately one hour. The grey-scale is the cross- 
power spectral density. The top row is the C0-C8 baseline and the 
bottom row the C0-W4 baseline, which are approximately 560 m 
and 9400 m respectively. Large bright patches indicate broadband 
RFI, and vertical lines indicate line RFI. The first panel in each 
row is the initial input data. The second is after the initial Scr 
mask, with most interference still visible. The third panel is after 
removing the largest eigenvalues in the SVD. The broadband in- 
terference is no longer visible. The last panel shows the final 3cr 
mask removing the line RFI, leaving a nearly uniform image. 



terference was removed. The first mask can not be too ag- 
gressive or the teclmiques to remove broadband interference, 
discussed below, will fail. 

Singular value decomposition (SVD) is used to separate 
broadband radio sources on tlie ground from those in the 
sky. Ground-based sources contribute most to the largest 
eigenvalues since they do not move as a function of time 
with respect to the array, while sky sources rotate. One 
iiour of data iias 14396 time records, and each record has 
about 7.5 million entries corresponding to the number of 
frequency cliannels and baselines between the 60 antennas. 
This is treated as a matrix with each time record as one 
row. The 50 largest eigenvalues are identified through an 
SVD and flagged as noise to be removed. A sample of the 
data at a few intermediate stages showing the successful re- 
moval of both line and broadband RFI can be seen in Fig. [2l 
The RFI patterns in (u, v) space both before and after the 
SVD are illustrated in Fig. [S] 



2.3 Physical Removal of Radio Interference 

Sources 

While some RFI can be removed in the data analysis, this 
risks removing sky signal as well. Ideally one would prevent 
RFI from occurring at all by identifying and correcting the 
physical sources. We can take advantage of the fact that as 
an interferometer GMRT is able to make images of the near- 




Figure 3. Raw visibilities with (u,v) distance |b| < 200 be- 
fore the SVD RFI removal step (left) and after the RFI removal 
(right). The bright patterns in the central vertical strip (low \u\) 
are caused by RFI, and are significantly reduced by the SVD pro- 
cedure. The effect of this process on the final power spectrum is 
considered in section [3.21 



field, and create maps of bright RFI sources near the GMRT 
itself. 

Candidate sources detected by a single baseline appear 
as a hyperbola of equal light arrival time in near-field images 
of a single SVD mode. When a source is detected by many 
baselines, the corresponding hyperbolas intersect at a single 
point in the image. These near-field images, an example of 
which is shown in Fig. [l] become the 'RFI maps' which 
are used to isolate interference sources. For sources with a 
high enough duty cycle, the GPS position given by the RFI 
maps are accurate enough that a handheld yagi antenna 
and portable radio receiver can be used to find the precise 
source. With the handheld antenna direction to a source 
can be determined by ear to better than 30 degrees, and 
thereafter localised by successive triangulation. 

By comparing the position of a noise transmitter rela- 
tive to the candidate interference source both physically and 
in the RFI map, we can confirm whether we have identified 
the correct source. Occasional misidentifications are likely 
due to the intermittent nature of the sources; there is no 
guarantee that the radiating source in the images will still 
be active when an attempt at identification is made. Since 
the beginning of this effort, calibration has been improved 
to locate sources within about 100 m. Sources could in prin- 
ciple be located using only the array with a precision 10 m, 
but the accuracy of available maps, GPS equipment, and 
other factors limit the real-world precision. 

Candidate RFI sources identified so far include trans- 
formers, power line junctions, and loose wires in contact with 
power lines. Table [1] lists some candidate sources recorded 
on 2008 December 6, and descriptions. In February 2009, 
GMRT began a collaboration effort with MSEDCL, which 
controls the power transmission lines in the area around 
GMRT, to assist us in removing local interference sources. 
A wire hanging over 132 kV liigh tension power lines (#1 in 
Table [l| which was by far the brightest interference source 
within 20 km, appearing in 30 of the 50 largest eigenmodes, 
was removed in February 2009. Other bright sources have 
also been removed since then. Although these sources are 
still present in the raw data from 2007, this procedure has 
shown that the SVD algorithm applied in the present work 
is successful in identifying real RFI sources during analysis. 
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Table 1. List of candidate RFI sources based on observations on 2008 December 6, identified by the strongest SVD mode number from 
which the GPS coordinates are derived. The two brightest sources, modes number 1 and 9, were removed in early 2009. 



Mode Longitude Latitude Description 



1 


74 


034653 


19 


147676 


Conclusively identified as a wire hanging over 132 kV high tension power lines, clearly visible 
using the handheld yagi antenna several kilometres away. Wire removed with assistance from 
MSEDCL on February 26, 2009. 


9 


74 


085335 


19 


171797 


Identified with wires hanging from an unused telephone pole, positioned directly under high 
tension lines. 


13 


74 


034271 


19 


037518 


Source in the area transmitted very intermittently, making identification difficult. Potentially 
a transformer 20 m north of coordinate. 


14 


74 


096581 


19 


209597 


Identified as a small wire hanging on a 500 kV power line, removed in February 2010. 


15 


74 


120018 


19 


181877 


Two possible sources, both transmission towers, separated by about 300 m. The one closest 
to the coordinate is a T-junction of two high voltage lines. 


17 


74 


072952 


19 


088457 


Initially identified to be a small pump approximately 150 m west of this coordinate, though 
other candidates include two nearby transformers. 


33 


74 


075623 


19 


108976 


Transformer approximately 100 m northwest, unambiguously radiating at low levels. 



2.4 Single-night Images 

After RFI removal has been completed in each one hour 
scan, polarization calibration is used to combine these into 
images of an entire night, typically of about 8 hours, as can 
be seen in Fig. [5] This step accounts for leakage between the 
left and right polarization signals, as well as effects caused 
by the relative rotation of the array and sky source by taking 
into account the changing parallactic angle with time. 

Although the array elements very nearly lie in a sin- 
gle plane, the rotation of the Earth turns this plane over 
time with respect to the line of sight. To correctly stack 
many hours together, one must consider the frequency and 
directional dependence of the measurement. One effect to 
consider is the change of the primary beam with frequency. 
For a given frequency, the primary beam is well defined and 
independent of (u, v). Deconvolution is difficult, but one can 
pick the frequency with the smallest beam and restrict the 
field of view at other frequencies to match. Computation- 
ally this can be done by convolving with the ratio of the 
beams, which is possible because of the commutativity of 
the primary beam operators. 

However, this is a small effect compared to the effects 
of a 10 term. Typically relative antenna positions are con- 
sidered to be on a two dimensional plane with coordinates 
parametrized by (w, v) in units of the observing wavelength. 
Since baselines become non-coplanar over the course of the 
night, a third dimension, w, must be included, without which 
the image becomes blurred. 

For a thorough exp l anatio n of w term issues, see for 



example ICornwell et al. (|2008|) a nd, for a more general 



overview, [Thompson et al] ( 200 ll ). As the observing fre- 
quency changes across the band, a single baseline samples 
different points in (u, v) space. At any given frequency, the 
(u, v) plane is typically sparsely sampled, but with many 
frequency channels there will often be data at a single (u, v) 
point for at least a few channels. One can subtract the aver- 
age from these, and under the assumption that foreground 
sources have the same spectral index as the calibrator, this 
removes the foregrounds very well. Unfortunately, when con- 
sidering the full (m, u, w) cube, the data is too sparse for this 
to work. Since the w coordinate is dependent on the position 
in (u, v) space, and can become very large, the effects of the 



w term can not be corrected at different frequencies in the 
same way as the beam. 

Strategies for correcting this effect for GMRT are be- 
ing developed, wh i ch build on those used in CMB studies 
iMvers et alj|2003l : iHobson fc Maisingeill2002l '). but we have 
not yet applied these to our data. In the interim we restrict 
ourselves to the short baselines for which w is small and such 
corrections are not needed. The short baselines are the ones 
most sensitive to the EoR signal, so we retain most of the 
EoR sensitivity of the array. 

The flux scale of the field is set by calibrating relative 
to the pulsar in the centre of the field for every baseline, 
receiver, and frequency. Since the pulsar flux is known to 
be variable, to set the absolute flux scale we identify the 
source with the highest fiux in the sky image, and set this 
peak value equal to the known flux of that source. The ex- 
act value found depends on the angular resolution, deter- 
mined by the maximum baseline length used. For all days 
analysed in this work, the brightest apparent source is the 
radio galaxy 5C 7.245. The fiux of this sou rce was meas ured 
as 215.1 ± 12.6 mJy at 1415 MHz by Willis et all (|l976 ') at 
the Westerbork Radio Telesco pe, and as 685 ± 47m Jy at 
408 MHz in the 5C catalogue jPearson fc Kuslll978l ') for a 
sp ectral inde x of 0. 93 ± 0.07. Following the same argument 
as iPen et al.l (|2009h . to achieve the best match with other 
bright sources in the field, we adopt a value of 1.6 Jy for 
5C 7.245. Although the fiux of this source changes by as 
much as 10 per cent across the band, this change is less 
than that due to the uncertainty in the spectral index. Ad- 
ditionally, tlu^_gaJajcylMstwo components separated by 
12arcsec (|Willott et al.ll200ll '). GMRT is capable of resolv- 
ing this structure with |b| > 2600, where b = {u,v) is the 
distance between the two antennas in the baseline in units of 
the observing wavelength. To calibrate the fiux we require a 
point source so that the fiux is not spread over multiple im- 
age pixels. This is achieved by limiting the maximum base- 
line length used while calibrating to |b| ^ 2000. At much 
smaller |b|, we become similarly limited by confusion. 

Removing the foregrounds ade quately is essentia l for 
detecting the EoR signal (see, e.g., IWang et al.l |2006| ) and 
much work has been d one on simul ating foregrounds (e.g. 
iBowman et al.i2 009: Jc lic et al.ll2010l) and designing removal 
strategies (|Morales et al.ll2006l : lHarker et aLll2009l : [Liu et all 
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Figure 4. (a) A near-field image, measuring about 40 km on each 
side and centred on GMRT. An RFI source identified through 
large SVD modes is clearly visible as a bright spot in the upper 
right, (b) A map of GMRT covering the same area as in (a). 
Each black circle indicates an antenna location, and the open 
circle corresponds to the RFI source, #14 in Table [T] identified 
in (a). 



|2009| . and others). For a review see lMorales fc Wvithel (|2010l ') 
and references therein. To remove foreground sources, we 
take advantage of the fact that the flux of such sources do 
not vary greatly with frequency. A signal originating from 
the EoR will appear as additional variation on top of the 
foreground signal. To model spectrally smooth sources, for 
every timestamp and baseline we subtract a piecewise linear 
fit between the median fluxes in frequency bins of 8, 2, or 




Figure 5. Dirty sky image made with 8 hours of data from 2007 
December 10. The maximum (u, v) distance is 2000 with grid size 
of approximately Ait = 8 giving a field of view of 7.3 degrees on 
each side. This same («, v) distance cut is used to calibrate the 
peak fiux of 5C 7.245 to 1.6 Jy for each day The RMS of this 
image is 44 m Jy within the primary beam. 



0.5 MHz. This is simple to implement and has the advan- 
tage over a polynomial fit of a having straightforward and 
local window function which is simple to interpret. Fore- 
grounds which do not vary over this range are removed, while 
features with variability at higher frequencies remain. This 
method results in an upper limit to the EoR signal since 
the measured power may still include residual foreground 
variation. By subtracting the mean flux over 2 MHz to re- 
move galactic foregrounds, noise levels can be lowered to 
about 2 mjy for maps of approximately 10 degrees. Bright 
point sources may play a significant ro le jPatta et al.ll201(]| : 
iLiu et al.ir2009l : lDi Matteo et al.ll2004r ). but are not treated 
separately here. 

To first order, this is the same as a boxcar average 
which we can write as 



(1) 



where D{iy) is the input data, D{u) is the data after filtering, 
and the window function is 



io(x) = 5{x) -H{x,Aiy). 



(2) 



Here, 5{x) is the Dirac delta function and H{x, Au) is a step 
function centred at x with width Au corresponding to our 
chosen filter, and of unit area. 

In Fourier space, this window function is 



sin(fc^Ai^/2) 
kuAu/2 ' 



(3) 



More precisely it is a boxcar average with a variable width, and 
therefore a somewhat different curvature. 



The 
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Figure 6. Fourier transform of our chiosen window function for 
Ai^ = 0.5 MHz (solid line). Dashed lines indicate where ui{fc||) = 
0.5, which sets the minimum fcy to which we are sensitive. This 
value scales inversely with Av. The peak of ?i)(A;|| )^/fc^ (the dot- 
dashed line) indicates the fcy to which we are most sensitive under 
the assumption that the power spectrum is proportional to . 

shown in Fig. [6l where we have denoted the wave number as 
kv to make exphcit that it is in units of inverse frequency. 
This can be rewritten in terms of k\\ using the fact that in 
our redshift range IMHz ~ 11.6h~^Mpc. When converted 
to units of /iMpc~^, we write the wave number as fey to 
emphasize that the window function acts on structure along 
the hne of sight only. 

After filtering in this way, we denote the for which 
w{ky') — 0.5 as corresponding to the minimum fcy along 
the line of sight to which we are sensitive, while smaller 
fcy are removed by the filter and will not contribute as 
strongly to the power spectrum. The three filters of 8 MHz, 
2 MHz, and 0.5 MHz correspond to minimum fcy of 0.04, 0.16, 
and 0.65/iMpc~^ respectively. Sin ce the power sp ectrum is 
thought to be proportional to k~^ jlliev et al.ll2008l '). the line 
indicates the fcy to which we are most sensitive. 

The 0.5 MHz foreground filter reduces the peak and 
RMS values by a factor of 50, as can be seen in Fig.[71 The fil- 
ter is most effective within the primary beam. Fig. [5]demon- 
strates the dominance of RFI in the filtered dirty maps, and 
the improvement that the SVD step provides. Without the 
SVD, the maps in Fig. [7| would have been a factor of 4 nois- 
ier. 



3 MULTIDAY ANALYSIS 

3.1 Differences Between Nights 

Since RFI and foreground signals are so much larger than the 
EoR signal, small errors in the subtraction of foregrounds 
could easily result in a spurious signal. To avoid this we use 
cross-correlations between multiple days to make a statisti- 
cal measurement. 

To gauge how successful a cross-correlation measure- 
ment might be, we are interested in the relative similarity 
of the different nights, which we can gauge by taking the 
difference of visibilities. Days which subtract well will show 
mostly noise, while days which subtract poorly will still show 
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Figure 7. Data from Dec 10. The top row is the sky image with 
a 11.4 degree field of view, with |b| < 200 and \u\ > 60 binned in 
the {u, v) plane by Au = 5. The bottom row is the visibilities in 
the same range with Au = 0.4 to show structure. The left column 
is before any foreground subtraction. In this image the dominant 
source is B2 0825-f 24 (or 4C 24. 17) just south of the FWHM of the 
primary beam, with a peak value of 2.2 Jy. RMS within the beam 
is 343 mjy. The right column is after a 0.5 MHz subtraction. The 
peak value of this image is 47 mJy with an RMS of 6.2 mJy, lower 
by a factor of about 50. If put on the same grey-scale as the image 
without foreground subtraction no features would be visible. The 
dominant source after this filter is 3C 200, well outside the beam. 




Figure 8. Dirty sky images before the SVD RFI removal step 
(left) and after the RFI removal (right) with |b| < 200, after the 
0.5 MHz foreground subtraction. The field of view is 11.4 degrees. 
RFI clearly dominates the map on the left, with peak value of 
73 mJy and an RMS of 18 mJy. After RFI removal, the peak drops 
to 32 mJy with an RMS of 4.2 mJy. 



evidence of bright sources. Visibilities are only used when 
data exists at the same position in both days. This serves as 
a test of the important of ionospheric fluctuations. Though 
the ionosphere is calibrated for along the hght of sight to the 
pulsar in the centre of the field, the ionosphere could change 
across the field of view. Pairs of days for which this change 
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Table 2. Peak flux and RMS of the difference of each two day pair 
available, in mjy, with a maximum (m, v) distance of 600, which 
corresponds to a maximum baseline of 1.2 km. All values are in 
mJy. To remove foregrounds, a 2 MHz linear filter was applied. 









Unflltered 


With 2 MHz 


fllter 


Subtracted pair 


Peak Flux 


RMS 


Peak Flux 


RMS 


Dec 10 


Dec 


11 


1856.1 


227.3 


137.9 


19.9 


Dec 10 


Dec 


14 


888.4 


185.6 


62.5 


11.6 


Dec 10 


Dec 


16 


278.4 


49.8 


27.1 


3.2 


Dec 10 


Dec 


17 


447.9 


64.9 


26.2 


3.6 


Dec 10 


Dec 


18 


548.3 


70.4 


32.5 


4.2 


Dec 11 


Dec 


14 


1037.8 


145.1 


256.1 


56.5 


Dec 11 


Dec 


16 


2148.3 


256.7 


158.1 


23.9 


Dec 11 


Dec 


17 


2528.2 


318.5 


106.1 


21.1 


Dec 11 


Dec 


18 


2997.5 


356.1 


159.1 


12.4 


Dec 14 


Dec 


16 


1221.5 


193.9 


50.2 


7.0 


Dec 14 


Dec 


17 


1311.1 


257.2 


50.4 


7.6 


Dec 14 


Dec 


18 


1433.5 


233.0 


67.1 


4.9 


Dec 16 


Dec 


17 


202.0 


78.9 


23.8 


2.8 


Dec 16 


Dec 


18 


224.2 


31.2 


24.4 


3.4 


Dec 17 


Dec 


18 


206.9 


60.6 


24.6 


3.2 



away from the field of view is different will not subtract well. 
The results of the subtractions are shown in Table (2] 

It can be seen in Table[2]that December 11 gives consis- 
tently poor results in both the unflltered and filtered images, 
and was thus excluded from all subsequent analysis. From 
these maps we can also conclude that RFI is not dominat- 
ing the differences from day to day. If this had been the 
case, since RFI is typically isolated in (u, v) space, it would 
be visible across the whole sky image, including far out- 
side the primary beam. However, it can be seen in Fig. [7] 
that structure decreases rapidly away from the centre of the 
field, indicating an astronomical source. This is true of all 
subtraction pairs. 

3.2 Cross-correlations and Power Spectra 

The power spectrum of sk y structure can be determ ined di- 
rectly from the visibilities (jZaldarriaga et aLll2004l ). To find 
the cross-power of two days, we take the product of the visi- 
bilities after gridding with a cell spacing equal to the size of 
the beam (Aw = 20). Calculating the power spectrum from 
the visibilities instead of from the correlations in the sky 
image takes advantage of the fact that the visib ilities have a 
nearly diagonal correlation matrix in the noise (i White et al.l 
I1999D . We then find the weighted average of visibilities in 
annuli of {u,v) space which gives the power in units of Jy^. 
Since the amount of data at large |b| decreases, we increase 
the width of each successive annuli by 60 per cent with in- 
creasing |b|. The smallest bin width is equal to the size of 
the beam. This prevents large artificial variability in power 
due to sparse sampling. Visibilities are weighted 'naturally', 
by the inverse of the noise. Since we expect our sensitivity 
to the EoR signal to diminish rapidly with increasing base- 
line length, we look only at the first few points, averaged 
over all possible cross-correlations. Additionally, it is known 
that the SVD will introduce a loss of power at low \u\. To 
avoid this, we impose a limit of \u\ > 60 when taking the 
cross-correlations, determined by requiring that the power 
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Figure 9. Comparison of the cross-power of December 10 with 
all other days under four different conditions. The dot-dashed 
and double-dot-dashed lines are before and after the SVD RFI 
removal, respectively, including all u. It can be seen that power is 
lost in the SVD. Similarly, the solid and dashed lines are before 
and after the SVD step, respectively, this time with a \u\ > 60 
limit imposed. In this case the two lines are almost identical, 
meaning the SVD had little effect on the total power. 

spectrum before and after the SVD differ by less than la 
as shown in Fig. |9l The part of the (u, v) plane that is lost 
with this cut can be seen in Fig. |3] 

The power spectrum of the cross-correlation can be con- 
verted to units of using 




where I = 27r|b|, b — {u,v) is the visibility coordinate in 
units of wavelength, Ci is the power measured i n K^, 6i, is 
the p rimary beam size, and u is the wavelength (|Pen et al.l 
[2003) ■ The quantity in the angle-brackets on the right is 
equal to the power in Jy^ found above. This conversion is 
written in terms of the GMRT observations with a primary 
beam of 6b — 3.3° and v — 150 MHz. If gridding in the 
(u, v) plane is too fine, the data become noisy, while very 
coarse gridding requires the assumption that the data are 
constant across the whole cell. As mentioned, we use a grid- 
ding equal to the size of the beam. A fully optimized es- 
timate would require a maximum likelihood code which is 
being worked on in the context of the Cosmic Background 
Imager (CBI) gridder ( Myers et al . 2003). Fig. [TUlshows the 
weighted-a verage of all cro s s-corr elation pairs, excluding De- 
cember ll. lBernardi et al] (|2009l ) reported a power spectrum 
of foregrounds without subtraction in the galactic plane at 
a level compara ble to our measurement both here and in 
iPen et al.l (|2009l ). 

W e have plott ed the power spectra with 2a bootstrap 
errors l|Efronlll979l ') which were derived as follows: Using five 
days of data, there are ten possible cross-correlations. From 
these, ten are randomly sampled, with replacement, result- 
ing in a slightly different power spectrum. This is repeated 
10* times, and the variance on this set of power spectra is 
calculated to give the error on the original. Formally this 
quantifies the error when taking independent samples of a 
statistical distribution. In our case this is not a rigorous er- 
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Figure 10. Average power spectrum in units of of all com- 
binations of days, excluding December 11, as a function of the 
multipole moment I. Each point is shown with a 2a upper limit 
derived from a bootstrap error analysis, which is in most cases 
smaller than the size of the point. The points are logarithmically 
spaced as described in the text, from left to right covering the 
ranges 377 < Z < 578, 578 < i < 899, and 899 <« < 1414. Tri- 
angles are the power before subtracting foregrounds, diamonds 
are after 8 MHz mean subtraction, squares are after 2 MHz mean 
subtraction, and circles are after 0.5 MHz s ubtraction. The c urved 
solid line is the theoretical EoR signal from ljelic et al.l l|2008l ) , and 
the dashed line is the theoretical EoR signal with a cold absorbing 
IGM as described in the text. 



ror, but a suitable straightforward estimate given the main 
complications discussed earlier. 



3.3 Comparison to Models 

Fig. [TO] can be compared to simulated result s from the 
Low F requency Array (LOFAR) EoR project in ljelic et al.l 
(|2008l '). which assumes Ts Tcmb- At low I, their simu- 
lated EoR signal is approximately (10 mK)^, while our low- 
est point with a similar 0.5 MHz bandwidth filter is (50 mK) 
with a 2a upper limit of (70 mK)^. These results are com- 
parable to the sensitivities LOFAR expects after 400 hours 
of the EoR project. We have also considered the case where 
reheating of the IGM does not occur, so the spin temper- 
ature renMjr]s_cougled^ kinetic temperature of the 
gas ijCiardi fc Madaul lioo^ t. In this case, the IGM cools 
adiabatically after decoupling from the CMB at z f» 150. 
The temperature fluctuations scale with (1 -I- Tcmb/Ts), 
and the power scales with the same factor squared. Using 
Tk = Tcmb(1 + z)/150 aX z = 8.6, the power becomes ap- 
proximately 275 times larger. This line is shown in Fig. 1101 
and is comparable t o the data. Scalin g up the warm IGM 
power spectrum from llliev" et al.1 (|2008l ) or ljelic et all (|2008l ') 
in this way is a reasonable approximation to what the ex- 
pected signal in a cold IGM. For a more detailed studies of 
th e signal in s u ch an abso rption regime, we point the reader 
to lBaek et all (|2009l ) and lBaek etaP (|2010l ). 

The pow er spectrum of reionizatio n is intrinsically three 
dimensional (|Morales fc HewittI |2004 ) . The strongest con- 
straints on the 3D power spectrum for the Az/ = 2 and 
0.5 MHz foreground filter case are shown in Fig. 1111 This 
uses A:^ = fen -f fc^, where fey is given by the windowing func- 



Figure 11. 3D power spectrum for the same data shown in 
Fig. 1101 using fc^ = fcy -|- fc^. This is dominated by fey, so the 
bin width in fej^ does not influence the horizontal position of the 
limits. The strongest constraints from the 2 MHz and 0.5 MHz fil- 
ters are shown (square and circle, respectively). Upper limits are 
2(T bootstrap errors. Three possible signals a re shown. The dashed 
line is the prediction from lliTev et al.ril2008l) and the double-dot- 
dashed line is the same for a cold IGM. The solid line comes from 
the single-scale bubble model as described in the text for a cold 
IGM, using fc = 2.5/-R to show the maximum power at all fc. For 
the two points shown, the bubble diametres which achieve this 
maximum power are 27 and 7.4/i~'^Mpc respectively. Only the 
0.5 MHz point imposes a limit on the diametre. For a warm IGM 
case, this signal would be reduced by the same factor as in the 
two dashed lines. 



tion of the filter and fcx = //6fe~^Gp c. Whe n comparing 
this to the prediction from llliev et al] lj2008h . one should 
note that our windowing function will also reduce the pre- 
dicted signal by at most a factor of two. 

We also consider an idealized case in which the ionized 
bubbles during reionization are of uniform scale and non- 
overlapping. Then for a given fc there will be a characteristic 
bubble radius R at which the power is maximized. By taking 
the 3D Fourier transform of a perfectly ionized bubble, and 
requiring that the universe is 50 per cent ionized, one can 
show the power is given by 



fc 



2-K 



cos fci? sin kR 



(kRY [kRf 



(4) 



and is maximized when kR ~ 2.5. In this case, 
fc^/(27r^)P(fc) -p^ T^/^, where Tb is the brightness tempera- 
ture ~ 30 mK in an X-ray heated IGM or almost -500 mK 
in a cold absorbing IGM. This signal would be more than 
an order o f mag nitu de larger t h an th e predictions by e.g. 
llliev et al.1 (120081 ') or Ijelic et all l|2008l ). In Fig. [11] we have 
included the power spectrum from this model with fc cho- 
sen to maximize the power in the range of interest. The 
data currently impose a limit on the size of bubbles in this 
single-scale model. Our upper limits with a 0.5 MHz fore- 
ground subtraction rule out bubbles with diametres from 
2.2 to 12.4/i~^Mpc in the redshift range 8.1 < z < 9.2. The 
cold IGM constraint is applicable even in the case of simula- 
tions, since only UV photons are included which themselves 
do not heat the IGM. 
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4 CONCLUSION 

The data analysis has been completed on six days from De- 
cember 2007 with a noise level of approximately 2 mjy on 
most nights. The SVD removal strategies for broadband RFI 
used lower noise by a factor of 4 in temperature, or 16 in 
power, which flagging alone cannot achieve. We have also 
tested for ionospheric variations and found that our pulsar 
calibration is sufficient for dealing with these effects. After 
RFI removal and foreground subtraction, we have measured 
a power spectrum which represents a new upper limit on the 
21 cm brightness temperature fluctuations during the epoch 
of reionization. These results can be used to constrain as- 
sumptions about the state of the IGM at these times, par- 
ticularly in the case of a Ly-a pumped, but cold, IGM. 

The previo us best limit o n 21c m signal at comparable 
redshift was by iBebbingtonI lll986h. who repo r ted n o fea- 
tures down to 5 K at 2 = 8.4. IParsons et aP (|2010( ) have 
reported a similar limit of about 5K using PAPER. The 
upper limit we present here is approximately 70 mK on the 
variance in 21 cm brightness temperature a.t z = 8.6, almost 
two orders of magnitude better than these previous limits. 
Residual foreground contamination and RFI may still be 
contributing to this power, but the EoR signal can not be 
larger than this. 

The analysis of additional observations of the B0823+26 
fleld made since 2007, containing approximately 2 times 
more data than has been treated in this work, is ongoing. 
We expect to continue to improve these results, and the 
planned addition of a second field will strengthen this sta- 
tistical measurement. The dominant uncertainty source re- 
mains RFI and foreground modelling. 

These results provide a flrst-look at the progress made 
at GMRT in detecting EoR. The GMRT-EoR team has been 
continuing observations, most recently with an additional 
150 hours allocated in the summer of 2010. We continue to 
improve both the system temperature of GMRT antennas 
and the RFI monitoring and removal strategies and expect 
to improve upon these results as the analysis of the newer 
data progresses. 



5 ACKNOWLEDGMENTS 

We thank Chris Hirata for his contributions to the analysis 
pipeline and a reading of the paper, and the anonymous re- 
viewer for many useful comments. We also thank the staff of 
the GMRT that made these observations possible. GMRT 
is run by the National Centre for Radio Astrophysics of 
the Tata Institute of Fundamental Research. The compu- 
tations were performed on CITA's Sunnyvale clusters which 
are funded by the Canada Foundation for Innovation, the 
Ontario Innovation Trust, and the Ontario Research Fund. 
The work of UP is supported by the National Science and 
Engineering Research Council of Canada. The map image 
in Fig. 2] is copyright by OpenStreetMajfl contributors and 
used under licensgj. 



http : //www ■ openstreetmap ■ org/] 

http : //creatlveconmions ■ org/llcenses/by-sa/2 ■ 0/ 1 



REFERENCES 

Alvarez M. A., Pen U., Chang T., 2010, ApJ, 723, L17 
Ananthakrishnan S., 1995, Journal of Astrophysics and As- 
tronomy Supplement, 16, 427 
Back S., Di Matteo P., Semelin B., Combes F., Revaz Y., 

2009, A&A, 495, 389 

Baek S., Semelin B., Di Matteo P., Revaz Y., Combes F., 

2010, preprint (astro-ph/1003.0834) 
Bebbington D. H. O., 1986, MNRAS, 218, 577 
Becker R. H. et al., 2001, AJ, 122, 2850 

Bernard! G., de Bruyn A. G., Brentjens M. A., Ciardi B., 
Harker G., Jelic V., Koopmans L. V. E., Labropoulos P., 
Offringa A., Pandey V. N., Schaye J., Thomas R. M., 
Yatawatta S., Zaroubi S., 2009, A&A, 500, 965 

Bowman J. D., Morales M. F., Hewitt J. N., 2006, ApJ, 
638, 20 

Bowman J. D., Morales M. F., Hewitt J. N., 2009, ApJ, 
695, 183 

Chen X., Miralda-Escude J., 2004, ApJ, 602, 1 
Chuzhoy L., Shapiro P. R., 2006, ApJ, 651, 1 
Ciardi B., Madau P., 2003, ApJ, 596, 1 
Cooray A., Li C, Melchiorri A., 2008, Phys. Rev. D, 77, 
103506 

CornweU T. J., Golap K., Bhatnagar S., 2008, IEEE Jour- 
nal of Selected Topics in Signal Processing, Vol. 2, Issue 
5, p.647-657, 2, 647 

Dalai N., Pen U., Seljak U., 2010, preprint (astro- 
ph/1009.4704) 

Datta A., Bowman J. D., Carilh C. L., 2010, ApJ, 724, 526 
Di Matteo T., Ciardi B., Miniati F., 2004, MNRAS, 355, 
1053 

Dijkstra M., Haiman Z., Loeb A., 2004, ApJ, 613, 646 
Djorgovski S. G., Castro S., Stern D., Mahabal A. A., 2001, 
ApJ, 560, L5 

Dunkley J., Komatsu E., Nolta M. R., Spergel D. N., Lar- 
son D., Hinshaw G., Page L., Bennett C. L., Gold B., 
Jarosik N., Weiland J. L., Halpern M., HiU R. S., Kogut 
A., Limon M., Meyer S. S., Tucker G. S., WoUack E., 
Wright E. L., 2009, ApJS, 180, 306 

Efron B., 1979, Ann. Statist., 7, 1 

Fan X., Narayanan V. K., Strauss M. A., White R. L., 
Becker R. H., Pentericci L., Rix H., 2002, AJ, 123, 1247 

Fan X., Strauss M. A., Becker R. H., White R. L., Gunn 
J. E., Knapp G. R., Richards G. T., Schneider D. P., 
Brinkmann J., Fukugita M., 2006, AJ, 132, 117 

Field G. B., 1959, ApJ, 129, 536 

Furlanetto S. R. et al, 2009, in astro2010: The Astronomy 
and Astrophysics Decadal Survey Vol. 2010 of Astronomy, 
Cosmology from the Highly-Redshifted 21 cm Line, pp 
82-+ 

Furlanetto S. R., Oh S. P., Briggs F. H., 2006, Phys. Rep., 
433, 181 

Furlanetto S. R., Zaldarriaga M., Hernquist L., 2004, ApJ, 
613, 1 

Gunn J. E., Peterson B. A., 1965, ApJ, 142, 1633 
Guo Q., Wu X., Xu H. G., Gu J. H., 2009, ApJ, 693, 1000 
Harker G., Zaroubi S., Bernardi G., Brentjens M. A., de 
Bruyn A. G., Ciardi B., Jelic V., Koopmans L. V. E., 
Labropoulos P., Mellema G., Offringa A., Pandey V. N., 
Pawlik A. H., Schaye J., Thomas R. M., Yatawatta S., 
2010, MNRAS, 405, 2492 



The GMRT-EoR Experiment: Hi Power Spectrum 11 



Harker G., Zaroubi S., Bcrnardi G., Brcntjens M. A., de 
Bruyn A. G., Ciardi B., Jclic V., Koopmans L. V. E., 
Labropoulos P., Mellema G., Offringa A., Pandey V. N., 
Schaye J., Thomas R. M., Yatawatta S., 2009, MNRAS, 
397, 1138 

Hobbs G., Lyne A. G., Kramer M., Martin C. E., Jordan 

C. 2004, MNRAS, 353, 1311 

Hobson M. R, Maisinger K., 2002, MNRAS, 334, 569 

Hogan C. J., Rees M. J., 1979, MNRAS, 188, 791 

Iliev I. T., Mellema G., Pen U., Mcrz H., Shapiro P. R., 

Alvarez M. A., 2006, MNRAS, 369, 1625 
Iliev I. T., Mellema G., Pen U.-L., Bond J. R., Shapiro 

P. R., 2008, MNRAS, 384, 863 
Jelic v., Zaroubi S., Labropoulos P., Bernardi G., de Bruyn 

A. G., Koopmans L. V. E., 2010, MNRAS, pp 1369-+ 
Jelic v., Zaroubi S., Labropoulos P., Thomas R. M., 

Bernardi G., Brentjens M. A., de Bruyn A. G., Ciardi 

B. , Marker G., Koopmans L. V. E., Pandey V. N., Schaye 
J., Yatawatta S., 2008, MNRAS, 389, 1319 

Kassim N. E., Lazio T. J. W., Ray P. S., Crane P. C, 

Hicks B. C, Stewart K. P., Cohen A. S., Lane W. M., 

2004, Planet. Space Sci., 52, 1343 
Komatsu E., Smith K. M., Dunkley J., et al 2010, preprint 

(astro-ph/1001.4538) 
Lidz A., Zahn O., McQuinn M., Zaldarriaga M., Hernquist 

L., 2008, ApJ, 680, 962 
Liu A., Tegmark M., Bowman J., Hewitt J., Zaldarriaga 

M., 2009, MNRAS, 398, 401 
Liu A., Tegmark M., Zaldarriaga M., 2009, MNRAS, 394, 

1575 

Loeb A., Zaldarriaga M., 2004, Phys. Rev. Lett., 92, 211301 
Lonsdale C. J. et al., 2009, IEEE Proceedings, 97, 1497 
Madau P., Meiksin A., Rces M. J., 1997, ApJ, 475, 429 
Mao Y., Tegmark M., McQuinn M., Zaldarriaga M., Zahn 

O., 2008, Phys. Rev. D, 78, 023529 
Masui K. W., McDonald P., Pen U., 2010, Phys. Rev. D, 

81, 103527 

McQuinn M., Lidz A., Zahn O., Dutta S., Hernquist L., 
Zaldarriaga M., 2007, MNRAS, 377, 1043 

McQuinn M., Zahn O., Zaldarriaga M., Hernquist L., 
Furlanetto S. R., 2006, ApJ, 653, 815 

Mezzacappa A., 2005, Annual Review of Nuclear and Par- 
ticle Science, 55, 467 

Morales M. F., 2005, ApJ, 619, 678 

Morales M. P., Bowman J. D., Cappallo R., Hewitt J. N., 
Lonsdale C. J., 2006, New Astronomy Review, 50, 173 

Morales M. F., Bowman J. D., Hewitt J. N., 2006, ApJ, 
648, 767 

Morales M. F., Hewitt J., 2004, ApJ, 615, 7 

Morales M. F., Wyithe J. S. B., 2010, ARA&A, 48, 127 

Myers S. T., Contaldi C. R., Bond J. R., Pen U., Pogosyan 

D. , Prunet S., Sievers J. L., Mason B. S., Pearson T. J., 
Readhead A. C. S., Shepherd M. C, 2003, ApJ, 591, 575 

Parsons A. R., Backer D. C, Foster G. S., Wright M. C. H., 
Bradley R. F., Gugliucci N. E., Parashare C. R., Bcnoit 

E. E., Aguirre J. E., Jacobs D. C, Carilli C. L., Herne D., 
Lynch M. J., Manley J. R., Werthimer D. J., 2010, AJ, 
139, 1468 

Pearson T. J., Kus A. J., 1978, MNRAS, 182, 273 

Pen U., Chang T., Hirata C. M., Peterson J. B., Roy J., 

Gupta Y., Odegova J., Sigurdson K., 2009, MNRAS, 399, 

181 



Peterson J., Pen U., Wu X., 2004, preprint (astro- 
ph/0404083) 

Pritchard J. R., Loeb A., 2010, Phys. Rev. D, 82, 023006 
Purcell E. M., Field G. B., 1956, ApJ, 124, 542 
Rottgering H. J. A. et al., 2006, preprint (astro- 
ph/0610596) 

Roy J., Gupta Y., Pen U., Peterson J. B., Kudale S., 
Kodilkar J., 2010, Experimental Astronomy, 28, 25 

Salvaterra R., Haardt F., Volonteri M., 2007, MNRAS, 374, 
761 

Scott D., Rees M. J., 1990, MNRAS, 247, 510 

Shapiro P. R., Ahn K., Alvarez M. A., Iliev I. T., Martel 

H., Ryu D., 2006, ApJ, 646, 681 
Sunyacv R. A., Zeldovich I. B., 1975, MNRAS, 171, 375 
Thompson A. R., Moran J. M., Swenson Jr. G. W., 2001, 

Interferometry and Synthesis in Radio Astronomy, 2nd 

Edition 

Trac H., Gnedin N. Y., 2009, preprint (astro-ph/0906.4348) 
Tseliakhovich D., Hirata C, 2010, preprint (astro- 
ph/1005.2416) 

Wang X., Tegmark M., Santos M. G., Knox L., 2006, ApJ, 
650, 529 

White M., Carlstrom J. E., Dragovan M., Holzapfel W. L., 
1999, ApJ, 514, 12 

Willis A. G., Oosterbaan C. E., de Ruiter H. R., 1976, 
A&AS, 25, 453 

Willott C. J., Delorme P., Omont A., Bergeron J., Delfosse 
X., ForvciUe T., Albert L., Reylc C, Hill G. J., Gully- 
Santiago M., Vinten P., Crampton D., Hutchings J. B., 
Schade D., Simard L., Sawicki M., Beelen A., Cox P., 2007, 
AJ, 134, 2435 

Willott C. J., Rawlings S., Blundell K. M., 2001, MNRAS, 

324, 1 

Wouthuysen S. A., 1952, AJ, 57, 31 

Zahn O., Lidz A., McQuinn M., Dutta S. a., 2007, ApJ, 

654. 12 

Zaldarriaga M., Furlanetto S. R., Hernquist L., 2004, ApJ, 
608, 622 

Zaroubi S., Silk J., 2005, MNRAS, 360, L64 



